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Fullerene peapods, that is carbon nanotubes encapsulating fullerene molecules, can offer enhanced 
functionality with respect to empty nanotubes. However, the present incomplete understanding of 
how a nanotube is affected by entrapped fullerenes is an obstacle for peapods to reach their full 
potential in nanoscale electronic applications. Here, we investigate the effect of Ceo fullerenes on 
electron transport via peapod quantum dots. Compared to empty nanotubes, we find an abnormal 
temperature dependence of Coulomb blockade oscillations, indicating the presence of a nanoelec- 
tromechanical coupling between electronic states of the nanotube and mechanical vibrations of the 
fullerenes. This provides a method to detect the Ceo presence and to probe the interplay between 
electrical and mechanical excitations in peapods, which thus emerge as a new class of nanoelec- 
tromechanical systems. 



Since their discovery 1 , carbon nanotubes encapsulat- 
ing arrays of Ceo fullerene molecules the so-called 
fullerene peapods, have attracted a lot of interest among 
physicists, chemists, as well as materials and electrical 
engineers. Peapods are considered as potential building 
blocks in nanoelectronics, offering enhanced functional- 
ity with respect to empty nanotubes. Their prospective 
applications include data storage devices^, single elec- 
tron transistors^"^, and spin-qubit arrays for quantum 
computing. Nevertheless, it is still unclear to what ex- 
tent the electronic properties of peapods, critical for their 
applications, are affected by the entrapped fullerenes. 

Theoretical studie d 10 " 12 ! predict dramatic changes in 
the electronic band structure of the nanotube, due to 
hybridization and corresponding shifts of SWNT en- 
ergy states and Ceo molecular orbitals. This is sup- 
ported by scanning tunneling microscopy experiments^ 
and photoluminescence excitation/emission mappin^^. 
On the other hand, no fundamental perturbation of the 
nanotube electron system has been found in photoe- 
mission experiments^ and electrical transport measure- 
ments on peapod quantum dots (QDs) 6 8 . In the latter 
case, clean electron transport with signatures of Kondo 
physics^ has been observed via single or double QDs, 
in the regime of strong electronic coupling V between 
energy states of the dot and its electrical leads. Co- 
herent transport has also been found for peapods filled 
with Sc@C82 metallofullerenes 17 , whose incorporation in 
principle should have an even stronger effect on the nan- 
otube's band structure than Ceo molecules. 

Here, we study the effect of the encapsulated Ceo 
fullerenes on electron transport via peapod quantum dots 
in the regime of weak electronic coupling. Compared to 
empty nanotubes, we find an abnormal temperature de- 
pendence of the Coulomb blockade (CB) conductance os- 



cillations, indicating a modification of single-electron ex- 
citations at low temperatures. These novel elementary 
excitations, which appear as a result of Ceo encapsula- 
tion, can be explained by a nanoelectromechanical cou- 
pling between quantized electronic states of the nanotube 
and the mechanical vibrations of the fullerenes. Their ap- 
pearance provides a method to detect the Ceo presence 
and to probe the interplay between electrical and me- 
chanical degrees of freedom in fullerene peapods, which 
emerge as a new class of nanoelectromechanical systems 
(NEMS). We stress that other types of devices for which 
a coupling to mechanical vibrations plays a crucial role 
in electrical transport are molecular conductors^ and 
transistors^, as well as suspended nanotube d 20 * 21 *. 

I. RESULTS 

Figure shows a transmission electron microscopy 
image of our peapod material, containing SWNTs 
densely filled with Ceo molecules. For electrical trans- 
port measurements, individual peapods with a diameter 
of 1.3-2 nm were contacted by Pd/Au source and drain 
electrodes, as schematically illustrated in Fig. [TJd. (See 
Methods for more details on devices.) The Ceo peapods 
showed clear Coulomb blockade behavior at low temper- 
atures (Fig. indicating that single-electron charging 
effects dominated the low-temperature electrical trans- 
port. Both the charging energy U ~ 8 — 12 meV and the 
electronic level spacing AE ^2 — 3 meV corresponded 
well to a 400 nm wide separation between the contact 
electrodes. Since the dots were probed at gate voltage 
ranges corresponding to a reduced electronic coupling T 
to the leads, no zero-bias Kondo anomaly^lwas observed 
here. 
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FIG. 1: Ceo fullerene peapod device, (a) Transmis- 
sion electron microscopy image of Ceo peapods with densely 
packed fullerene chains. The scale bar is 5 nm long, (b) 
Schematic layout of the device. An individual Ceo peapod de- 
posited on a Si02 substrate is contacted by Pd/Au source and 
drain electrodes, spaced 400 nm apart, (c) Stability diagram 
of a Ceo peapod. The (color coded) differential conductance 
dl/dV is measured at T = 0.3 K as a function of backgate 
voltage V g and bias voltage V, showing the diamond-shaped 
features characteristic for the Coulomb blockade. 



Figure [2^i depicts the linear conductance G of a 
fullerene peapod as a function of gate voltage V g , reveal- 
ing Coulomb blockade oscillations. This kind of experi- 
ment is a basic tool for probing both the energy spectrum 
and the wavefunctions of the elementary excitations in a 
quantum dot: The positions of the CB peaks provide in- 
formation about the spectrum, while the magnitude of 
the tunneling current is proportional to the overlap be- 
tween the wave functions of the incident electron and the 
elementary excitation in the dot. The temperature de- 
pendence of the heights of the five CB peaks that appear 
in Fig. [2^i is shown in Fig. [2}d. 

The lineshape and temperature response of Coulomb 
blockade oscillations in the linear conductance has been 
widely investigated in the past!^, both for metallic and 
semiconductor quantum dots. In the CB regime where 
electron transport occurs only via a single electronic dot 
level, the tunneling probability as a function of energy is 
given by the Breit-Wigner formula (see caption to Fig.j2|. 
The temperature dependence of the conductance peaks 
is therefore entirely due to thermal broadening of the 
energy of the incident electrons: The peak height is in- 
dependent of temperature in the low-temperature limit, 



where k^T <C T, while in the opposite high-temperat 
limit, k^T ^> T, a universal 1/T-behavior appears. This 
corresponds to what we observed for empty SWNTs, as 
shown in Figs. [2]>e. 

For fullerene peapods, however, the measured temper- 
ature dependence of the conductance peaks is often sig- 
nificantly different (see Fig. 2b), implying that the Breit- 
Wigner picture is no longer valid. We interpret this dis- 
crepancy as a signature of a change in the spectrum of 
the elementary excitations and of their wavefunctions, 
a change that can be associated with the encapsulated 
fullerenes. Typically, as highlighted in Fig.[3^i, an anoma- 
lous temperature dependence appears on the scale of a 
few kelvin. The corresponding energy is of the order of 
0.1 meV, which is small on the electronic energy scale but 
relevant for the longitudinal vibrational motion of the 
encapsulated fullerenes. (For the low-diameter peapods 
used in our study, the characteristic scale of transverse 
vibrations would be ~1 meV.) 

In order to interpret our experimental results, we use 
a model where SWNT electrons are electromechanically 
coupled to a set of low- frequency vibrators: single Ceo 
fullerene^ or their clusters trapped by a longitudinal 
harmonic potential inside the nanotube. Within this ap- 
proach, the peapod is treated as a multilevel quantum 
dot for which each quantized electronic level is selectively 
coupled to a single vibrational mode, in a similar manner 
as in the local Holstein modePS, and the corresponding 
dimensionless electron- vibron coupling constant A is pro- 
portional to the amplitude of the vibronic quantum fluc- 
tuations. As a result, the resonant tunneling electrons 
can be dressed by vibrational excitations if they spend 
sufficiently long time in the dot to form real electron- 
vibron hybrids - Holstein polarons^. Since the proba- 
bility of polaron formation is enhanced for long electron 
dwell times in the quantum dot, 1/T, polaronic effects be- 
come more pronounced in the regime of weak electronic 
coupling T. 

We stress that the electron-vibron coupling does not 
have any significant effect on the fullerene vibrations 
themselves, as they are mostly determined by the Ceo 
confinement within the carbon nanotube. Even though 
the coupling to each single electronic level may be strong 
and depend on the position of the oscillating fullerenes, 
their total energy shift averages out over a large number 
of individual electronic states in the dot. 

Our model calculations (Fig. ^$jp) have revealed quali- 
tatively similar features to those experimentally observed 
for conductance peaks deviating from conventional Breit- 
Wigner behavior (Fig. [3^i). In particular, three charac- 
teristic thermal regimes can be distinguished: (i) When 
the thermal energy is much smaller than the vibration en- 
ergy, k^T <C Tiujq, the conductance becomes suppressed 
for strong electron-vibron coupling A > 1. Such a po- 
laronic (or Franck- Condon) blockade of conductance^ 
signals the reduced probability for an electron to tunnel 
from a pure electronic state in the lead to a polaronic 
state in the dot, and vice versa, (ii) As the temperature 
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FIG. 2: Temperature dependence of linear conductance, (a) Coulomb blockade oscillations in the linear conductance 
G = dI/dV\v=o, measured as a function of gate voltage V g for a Ceo fullerene peapod. (b) Maximal conductance amplitude 
Gmax vs temperature T, determined for the five peaks indicated in (a). Peaks P4 and P5 behave as predicted by the Breit- 
Wigner (BW) model, while the others show an abnormal temperature dependence, (c) and (d) Lineshape of a conductance peak 
for a reference SWNT device, measured at (c) 0.3 K, and (d) 9.1 K. These lineshapes correspond to a thermally broadened 
Breit-Wigner peak G(E) = f dE'(-df(E')/dE')G B w(E - E'). At low temperatures, k B T < T, the peak approaches the 
Lorentzian lineshape Gbw(^) oc 1/[(E — Eq) 2 + (r/2) 2 ], as becomes apparent in (c). At higher temperatures, Y <C k B T, it can 
be approximated by a negative derivative of the Fermi function, G(E) oc —df/dE, see (d). (e) C ma x versus T for the reference 
SWNT device. As predicted by the BW theory, G max (T) is nearly independent of T at low temperatures, k B T <C T, while 
at higher temperatures, k B T <C T, the 1/T behavior is clearly recovered. In (b)-(e), the thin solid lines indicate fits to the 
experimental data performed using the Breit-Wigner model. 



rises, the population of the thermally excited vibronic 
levels increases which has a destructive effect on the po- 
laronic wavefunction. This leads to a non-monotonic be- 
havior of the conductance in the temperature interval 
/cbT < \ 2 huJo, where X 2 huJo represents the character- 
istic energy of a polaron. (iii) At high temperatures, 
/cbT ^> A 2 fk^o, polaronic effects are lifted and the con- 
ductance takes standard Breit-Wigner form G max oc 1/T. 

To determine the electron- vibron coupling A in the pea- 
pod devices, we have numerically fitted the temperature 
dependence of conductance maxima (Fig. [3^t) with the 
thermal progression predicted by our model (Eq. [5] in 
Methods). The obtained fits are in good agreement with 
the experimental data, yielding large values for the po- 
laronic coupling parameter, A = 1.6 — 1.9. These values 
are consistent with the suppression of the conductance 
observed at low temperatures. They are also of the same 



order of magnitude as in previ ous ex periments on sus- 
pended carbon nanotube devices^ 2 ^ where the vibronic 
states were due to intrinsic modes of the SWNT. (Note 
that in Refs . l20l and [2T1 the coupling constant was defined 
as g = A 2 ). The inset to Fig. |3^i gathers all the fitting 
parameters extracted for the investigated peaks. 

The vibrational origin of the abnormal features ob- 
served in the G max (T) is in addition supported by the 
appearance of weak, quasi-periodic excitation lines along 
the edges of the Coulomb blockade diamonds in the sta- 
bility diagrams, as indicated by arrows in Fig. The 
mean energy spacing ~ 0.25 meV of these nearly equidis- 
tant excitations is an order of magnitude smaller than the 
electronic level spacing AE = 2 — 3 meV in our 400 nm 
long peapod quantum dots. On the other hand, it closely 
matches the values of vibrational quanta extracted from 
the fits to G max (T), see the inset to Fig. [3^1. 
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FIG. 3: Polaronic interactions in a fullerene peapod. (a) G max vs T for conductance peaks with an abnormal temperature 
behavior with respect to the Breit-Wigner model. Solid and open symbols indicate the data points measured for two different 
peapod devices. The solid lines are fits to the experimental data, calculated using the polaronic model described in the text. The 
inset gathers the determined values of the fitting parameters, (b, Inset) A schematic diagram of the vibronic confining potential. 
The black parabola represents the harmonic potential defining the oscillating motion of the vibrator, while horizontal lines are 
vibronic energy states spaced by Tiujq. At temperatures k&T <C Tiujq, only the vibronic ground state (solid line) participates in 
the conductance. As soon as k^,T becomes of the order of Tiojq, thermally excited vibronic modes of the oscillating scatterer 
are generated (dashed lines) with characteristic energy n(z)hu)o, where n(z) is the Bose distribution with z = huo/k-BT. The 
ratio between this energy and the polaronic energy shift \ 2 Tiujq determines how detrimental the thermal vibrations can be to 
the polaronic effects, (b, Main figure) Calculated renormalized peak conductance G ma x = Gmax^o/7 as a function of reduced 
temperature T = k^T /Tiujq for three different values of the polaronic coupling constant A. At low temperatures, k^T <C Tkjo, 
an exponential suppression of the conductance occurs, G ma x oc G m &x,\=o(T) exp(— A 2 ), due to a reduced overlap between the 
electronic wave function in the lead and the polaronic wave function in the dot. Here, G ma x,\=o(T) is the thermally broadened 
BW peak conductance for noninteracting electrons. The stronger the polaronic coupling A, the stronger the suppression of the 
conductance is. By increasing the temperature, the competition between the thermally generated vibrons and the polarons 
gives rise to a non-monotonic behavior. An abnormal temperature dependence of the peak height is particularly well visible 
for /cbT < X 2 huJo where the conductance is dominated by the term (1/a/T) exp(— \ 2 Tiujq/ Ak^T). At higher temperatures, 
/cbT ^> A 2 /kJo, the polaronic blockade is completely overcome by thermal excitations and the conductance scales as 1/T, 
following the Breit-Wigner formula. See Methods section for more details on the model, (c) Closer look at a stability diagram 
dl/dV vs V vs V g . Arrows indicate weak quasi-periodic excitation lines running parallel to the Coulomb diamond edge. The 
inset shows their excitation energies; the solid line has a slope of 0.25 meV. 



II. DISCUSSION abnormal features in their temperature dependence, see 

Fig.|2]3. This observation can be explained if, for different 
SWNT electron states n, there are strong fluctuations^ 
As mentioned before, our experimental data clearly in- in the value of the polaronic coupling constant. Such fluc- 
dicate that not all of the studied conductance peaks show tuations may be due to a local character of the electrome- 
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chanical coupling. Consequently, only for some electronic 
levels the characteristic time needed to form a polaronic 
state l/(X^ l huJo) is sufficiently shorter than the electron 
dwell time in the dot 1/T n . 

We note here that the electronic coupling T could 
not be directly determined from the polaronic fits. In- 
stead, we were able to probe a related parameter 7 = 
r L lV(r L + T R ), yielding 7 = 0.004 - 0.050 meV for 
the investigated conductance peaks. Here, Tl and Tr 
are electronic couplings to the left and to the right reser- 
voir of the quantum dot, respectively, and r = Tl + Tr. 
Assuming symmetric barriers, Tl = Tr, the value of T 
could thus be evaluated as T = 47. For the asymmetric 
coupling, = 10r#, this would yield T = 12. 17. 

In order to gain some insight into the nature of the 
vibrational motion of fullerene molecules trapped inside 
a SWNT, we have performed a simple estimation of the 
potential that confines their longitudinal motion. The 
confining potential can readily be obtained from Giri- 
falco's universal potential for graphitic structures^ if we 
assume that it is mainly due to the interaction of a vi- 
brating Ceo molecule with its nearest-neighbor molecules 
(one on each side). The resulting potential is well ap- 
proximated by a harmonic potential characterized by the 
force constant k = 2.2 Nm" 1 . The corresponding vibra- 
tional energy quantum Tiujq matches the value of about 
0.25 meV, needed to fit the experiments, for a particle 
mass about ten times larger than that of a single Ceo 
molecule (Mc 60 = 1.2 • 10 -24 kg). The implication is 
that the relevant molecular vibrations may involve clus- 
ters of about ten Ceos. For peapods with densely packed 
fullerene chains, clusters of this size could form inside the 
SWNT, for example, as a result of local tube corrugations 
or presence of defects (openings) on SWNT side walP. 

The present incomplete understanding of how a nan- 
otube is affected by encapsulated fullerenes is an obstacle 
for peapods to reach their full potential in nanoelectronic 
applications^^. The experiments presented here shed 
light on this issue: They demonstrate the presence of na- 
noelectromechanical coupling in carbon peapods which 
thus emerge as a new class of nanoelectromechanical 
systems. Moreover, the characterization technique de- 
scribed here can be thought of as a tool for detecting 
the presence of fullerenes, or other molecules, in a nan- 
otube. Such an experimental verification can be per- 
formed in situ for individual tubes in electronic devices, 
unlike in other detection techniques which require either 
a large ensemble of peapods (Raman spectroscopy^, x- 
ray^) or probing in dedicated systems (transmission elec- 
tron microscopy, scanning tunneling microscopy 13 , low- 
temperature atomic force microscop}^) . 

Further studies using suspended peapods could allow 
a tuning of the fullerenes' spatial distribution within 
the nanotube by means of mechanical bending deforma- 
tions of the SWNT. If the resonant character of electron 
transmission could be controlled mechanically by these 
means or others, carbon peapods could find applications 
as NEMS resonators with an operation frequency set by 



the vibrational motion of the encapsulated fullerenes in 
the 100 GHz range. 



III. METHODS 
A. Devices 

The high-quality fullerene peapods investigated in our 
study were based on arc discharge-grown SWNTs. The 
empty tubes (also used as the reference SWNTs) were 
supplied by NANOCARBLAB as a purified material 
(80% pure grade). Due to the purification procedure 
(multi-step heating in air, soaking in HNO3, microfil- 
tration), the as-received SWNTs were already opened 
and required no further treatment prior to their filling 
with C 60 fullerenes (98% purity, INTERCHIM). The fill- 
ing was carried out via the vapor phase methocP, in an 
evacuated and sealed ampoule at 500° C for 24 hours. 
In order to remove excess fullerenes, heat treatment was 
afterwards performed under dynamic vacuum, at 800° C 
for 1 hour. Transmission electron microscopy indicated 
that 90-95% of the tubes in the resulting material en- 
capsulated fullerenes (i.e., were peapods). A quantita- 
tive analysis by x-ray diffraction^ confirmed that the 
peapods were densely packed with Ceos, with a filling ra- 
tio of ~ 80%. In addition, signatures of the encapsulated 
Ceo fullerenes were found in Raman spectroscopy. 

The peapods (or the reference empty nanotubes) were 
sonicated in dichloroethane and deposited on a silicone 
wafer covered with silicone oxide. Individual tubes with 
a diameter of 1.3 — 2nm were located and selected for 
further processing by means of atomic force microscopy. 
Electron beam lithography was then used to define source 
and drain electrodes of Pd/Au (15nm/45nm), with a 
spacing of 400 nm. A highly n-doped Si substrate, sit- 
uated beneath the 500nm-thick Si02 layer, served as 
a backgate. Electrical transport experiments were per- 
formed in a 3 He cryostat with a base temperature of 
T = 300 mK. 



B. Model 

We highlight here the main steps of our theoretical 
modeling, which is based on the local Holstein model 24 . 
More details can be found in Ref. [23] . 

In order to get analytical results, we assume that only 
one vibrating scatterer (a cluster of fullerenes or a single 
molecule) is confined inside the SWNT. The full Hamil- 
tonian for a two-terminal peapod QD system reads: 

T~L = i^L + #R + -f^T,L + #T,R + ^QD, (1) 

where H a = £k a c k c k a is the Hamiltonian for the 
Fermi leads and a = L, R for the left and right reservoir, 
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respectively. The peapod-QD Hamiltonian reads: 



where 



QD 



E 



V n (L,l)cic n (tf + b) + hw tfb, 

(2) 

where {s n } is the set of energy levels in the dot, 
hojQ is the vibrational quantum of energy, V n (L,l) = 
hvFXo(dk n /dx), and xq is the amplitude of the zero-point 
fluctuations of the bosonic field. c n (c}J and bQft) are 
fermionic and bosonic operators with canonical commu- 
tation relations. The longitudinal motion of the scatterer 
in the SWNT is defined by a harmonic potential. Finally, 
the tunneling Hamiltonian reads: 



H T ,a = ^t(4 a c n + h.c.) 



(3) 



where we suppose that the hopping matrix elements t 
between the two leads and the dot are of the same order. 

The above QD Hamiltonian can be diagonalized by 
a unitary transformation, which reveals that the major 
effect of the vibrational bosonic field is to renormalize 
the energy levels in the dot. These levels are shifted by 
an amount A^TkJo, the so called polaronic shift, where 
A n = V n (L,l)/(huJo) is the electron-vibron coupling con- 
stant. To obtain the expression for the conductance of 
our system, we start from the Meir-Wingreen formula^ 



I(V) 



2e 



de [h(e) - / R (e)]3m{Tr[ 7 G$ ) ( e )]}, 



(4) 

where /l/r(£) are the Fermi functions for the left and 
right reservoirs, 7 = r L r R /(r L + T R ); r L/R = 27rt 2 p L/R 
are the widths of the resonance due to tunneling from the 
left or right lead, and Pl/r are the densities of states in 
the reservoirs. Gqd(£) is the retarded Green function of 
the dot which can be evaluated in the wide-band limftP^. 

By means of the gate voltage, an energy level in the 
QD can be tuned to line up with the Fermi level £p in 
the leads, so that Sp = e n — X^Tiuoq. This is when the 
conductance reaches its maximum. In the linear conduc- 
tance regime, one finds that the temperature dependence 
of these maxima is given by the simple expression 



G max ,A(D = C maX) A=o(^)^A(^), 



(5) 



G 



max, A: 



-o(T) 
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2 h k B T 



(6) 



is the standard resonance conductance of a single-level 
quantum dot. In equation ([5|, the function F\(z), with 
z = hoJo/k-BT, can be expressed as a series, 



Fx(z) 



X 



exp{-A 2 [l 
exp(- 



f 2n(z)}} 

■^/2)/ ; [2A 2 v / ^)(l~ 



n 



(*))] 



(7) 

cosh^z/2) 

where 1\ is the modified Bessel function of the first kind 
and n(z) is the Bose distribution function. The thin solid 
lines in Fig. [3^i show fits to the experimental data ob- 
tained with the help of expression ([5| for the conduc- 
tance. 
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